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We report on a measurement of the inclusive jet production cross section in pp collisions at 
y/s = 1.96 TeV using data collected with the upgraded Collider Detector at Fermilab in Run II 
(CDF II) corresponding to an integrated luminosity of 385 pb” 1 . Jets are reconstructed using the 
kx algorithm. The measurement is carried out for jets with rapidity 0.1 < |y jet | < 0.7 and trans¬ 
verse momentum in the range 54 < p^ 4 < 700 GeV/c. The measured cross section is in good agree¬ 
ment with next-to-leading order perturbative QCD predictions after the necessary non-perturbative 
parton-to-hadron corrections are included. 

PACS numbers: PACS numbers 12.38.Aw, 13.85.-t, 13.87.-a 


The measurement of inclusive jet production in pp col¬ 
lisions at y/s = 1.96 TeV constitutes a test of perturba¬ 
tive QCD (pQCD) t l] predictions over more than eight 
orders of magnitude in cross section, and is sensitive to 
the presence of physics beyond the standard model [2j. 


The increased center-of-mass energy and integrated lu¬ 
minosity in Run II at the Tevatron make it possible to 
measure the cross section for jets with transverse mo¬ 
mentum p], p J T ct , up to about 700 GeV/c, extending the 
p^ 4 range by more than 150 GeV/c compared with pre- 
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vious results Jj]. This letter presents a new measure¬ 
ment of the inclusive jet production cross section as a 
function of p ^ 4 for jets with p ^ 4 > 54 GeV/c and ra¬ 
pidity [ 5 ] in the region 0.1 < |yJ et | < 0.7, where jets 
are reconstructed with the kx algorithm Hi Similar 
measurements have been carried out using cone-based jet 
algorithms in Run II Jjj. However, the kx algorithm has 
been widely used for precise QCD measurements at both 
e + e _ and e ± p colliders, and allows a well defined compar¬ 
ison to the theoretical predictions, without introducing 
additional parameters [ 6 ] in the calculations. The mea¬ 
surements are corrected for detector effects back to the 
particle (hadron) level and compared to pQCD next- 
to-leading order (NLO) predictions p|. A proper compar¬ 
ison with the theory requires corrections to account for 
non-perturbative contributions that become important 
at low pi/ 4 . This could explain the marginal agreement 
observed in previous studies TO]. 

The CDF II detector is described in detail else¬ 
where m. Here, the sub-detectors most relevant for 
this analysis are briefly discussed. The detector has a 
charged particle tracking system immersed in a 1.4 T 
magnetic field, aligned coaxially with the beam line. A 
silicon microstrip detector provides tracking over the ra¬ 
dial range 1.35 to 28 cm and covers the pseudorapid¬ 
ity 1 range \tj\ < 2. A cylindrical 3.1 m long open-cell 
drift chamber covers the radial range from 44 to 132 cm 
and provides full tracking coverage for \rj\ < 1. Seg¬ 
mented sampling calorimeters, arranged in a projective 
tower geometry, surround the tracking system and mea¬ 
sure the energy flow of interact ing particles in |? 7 | < 3.6. 
The central barrel calorimeter jl2|, [l^J covers the region 
1 77 1 < 1. It consists of electromagnetic and hadronic 
calorimeters segmented into 480 towers of size 0.1 in 
77 and 15° in <fi. The measured energy resolution for 
electrons is © 2%. The single-pion 

energy resolution, as determined from test-beam data, 

is . 50% © 3%. A hadronic calorimeter comple- 

y/E T [GeV] 

ments the coverage of the central barrel calorimeter in 
the region 0.6 < \r)\ < 1.0 and provides additional for¬ 
ward coverage out to \r]\ < 1.3. The forward region, 
1.1 < \r]\ < 3.6, is covered by scintillator-plate electro¬ 
magnetic and hadronic calorimeters. Cherenkov counters 
in the region 3.7 < \r]\ < 4.7 measure the number of in¬ 
elastic pp collisions to compute the luminosity | 3 | ■ 

Monte Carlo event samples are used to determine the 
response of the detector and the correction factors to the 
hadron level. The generated samples are passed through 
a full CDF detector simulation (based on GEANT3 |n| 
where the gflash Jn| package is used to simulate the 
energy deposition in the calorimeters) and then recon¬ 
structed and analyzed using the same analysis chain as 
for the data. Samples of simulated inclusive jet events 
have been generated using the PYTHIA 6.203 na and 
herwig 6.4 [i]| Monte Carlo generators. CTEQ5L ilg] 


parton distribution functions (PDFs) are used for the 
proton and antiproton. The PYTHIA samples have been 
created using a special tuned set of parameters, denoted 
as PYTHIA-TUNE A, that includes enhanced contributions 
from initial-state gluon radiation and secondary parton 
interactions between proton/antiproton beam remnants. 
Tune A was developed with dedicated studies of the un¬ 
derlying event using the CDF Run I data 2Q], and de¬ 
scribes the measured jet shapes in Run II |21(. 

The kx algorithm is used to reconstruct jets from the 
energy depositions in the calorimeter towers with trans¬ 
verse momentum above 0.1 GeV/c. First, all towers are 
considered as protojets. The quantities kx,i = Px; and 
kp,(i,j) = min(px ; , Pxj)-AR?j/D 2 , are computed for each 
protojet and pair of protojets respectively, where px,i de¬ 
notes the transverse momentum of the i th protojet, A 
is the distance (y — <f> space) between each pair of proto¬ 
jets, and D is a parameter that approximately controls 
the size of the jet. All kx,i and k T ,(i,j) values are then 
collected into a single sorted list. In this combined sorted 
list, if the smallest quantity is of the type kx,i, the cor¬ 
responding protojet is promoted to be a jet and removed 
from the list. Otherwise, if the smallest quantity is of 
the type kx.pj), the protojets are combined into a single 
protojet by summing up their four-vector components. 
The procedure is iterated over protojets until the list is 
empty. The jet transverse momentum, rapidity, and az¬ 
imuthal angle are denoted as PxVjal 1 Ycal’ an< ^ ^CAm 
respectively. The same jet algorithm is applied to all 
the final-state particles in the Monte Carlo to search for 
jets at the hadron level. The resulting hadron-level jet 
variables are denoted as Px'hadj ^had> an d '/'had- 

The measurements presented in this letter correspond 
to a total integrated luminosity of 385 ± 22 pb” 1 . Events 
are selected online using three-level trigger paths 0 
with different prescales. In the first-level trigger, a sin¬ 
gle trigger tower with Et above 5 GeV or 10 GeV is 
required. In the second-level trigger, clusters are formed 
around the selected trigger towers, and a cluster with Et 
above 15 to 90 GeV, depending on the trigger path, is re¬ 
quired. In the third-level trigger, jets are reconstructed 
using a cone-based algorithm, and a jet with Et above 20 
to 100 GeV is required. Jets are then searched for using 
the kx algorithm, as explained above, with D = 0.7. For 
each trigger data sample, the threshold on the minimum 
Pt'cal i s chosen in such a way that the trigger selection 
is fully efficient. The events are required to have at least 
one jet with rapidity in the region 0.1 < |ycAiJ < 0.7 
and corrected transverse momentum (see below) above 
54 GeV/c. The events are selected to have at least one re¬ 
constructed primary vertex with z-position within 60 cm 
around the nominal interaction point. In order to re¬ 
move beam-related backgrounds and cosmics rays, the 
events are required to fulfill I^t/V^Et < -F(Px^GiLL jet ), 
where 1 /t denotes the missing transverse energy 12 .‘jt | and 
YiEt = Y2i Et the total transverse energy of the event, 
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as measured using calorimeter towers with E l T above 
100 MeV. The threshold function -^(PxJcal je *) defined 

as F(p J ^) = min(2 + 0.0125 x pi^, 7), where PxJcal jet 
the uncorrected transverse momentum of the leading jet 
(highest p-J^) and the units are GeV. This criterion is de¬ 
signed to preserve more than 95% of the QCD events, as 
determined from Monte Carlo studies. A visual scan for 
P't'cal > 400 GeV/c showed no remaining backgrounds. 

The measured jet transverse momentum includes addi¬ 
tional contributions from multiple proton-antiproton in¬ 
teractions per bunch crossing at high instantaneous lu¬ 
minosity. The data were collected at instantaneous lu¬ 
minosities varying between 0.2 x 10 31 cm~ 2 s _1 and 9.6 x 
10 31 cm _ 2 s _1 . At the highest instantaneous luminosities, 
an average of two interactions per bunch crossing are 
produced. This mainly affects the measured cross sec¬ 
tion at low Px* where the contributions are sizeable. In 
CDF, multiple interactions are identified via the pres¬ 
ence of additional primary vertices reconstructed from 
charged particles. The measured jet transverse momenta 
are corrected for this effect by removing a certain amount 
of transverse momentum, e, for each additional primary 
vertex. A value e = 1. 62+^40 GeV/c is determined from 
the data by requiring that, after the correction is applied, 
the ratio of cross sections at low and high instantaneous 
luminosities does not show any p^* dependence. 

The reconstruction of the jet variables in the calorime¬ 
ter is studied using Monte Carlo. These studies indi¬ 
cate that the angular variables of a jet are reconstructed 
with no significant systematic shift and with a resolu¬ 
tion better than 0.05 units in y and <j> at low P'x'cal : 

improving as P't'cal increases. The measured jet trans¬ 
verse momentum systematically underestimates that of 
the hadron level jet, which is mainly attributed to the 
non-compensating nature of the calorimeter |24|. For jets 
with Px*cal nbout 50 GeV/c, the jet transverse momen¬ 
tum is reconstructed with an average shift of —19% and 
a resolution of 14%. The jet reconstruction improves as 
P't'cal increases. For jets with PxVjal about 500 GeV/c, 
the average shift is —5% and the resolution is about 7%. 
The bisector method j2f| is employed to evaluate how 
well the Monte Carlo reproduces the measured jet en¬ 
ergy resolutions. Data and Monte Carlo agree within a 
relative uncertainty of ± 8 % over the whole P'xVjal range. 

The measured PxVjal distribution is unfolded back 
to the hadron level using Monte Carlo event samples. 
PYTHIA-TUNE a provides a reasonable description of the 
different jet and underlying event quantities, and is used 
to determine the correction factors in the unfolding pro¬ 
cedure. In order to avoid any bias on the correction fac¬ 
tors due to the particular PDF set used, which trans¬ 
lates into slightly different simulated Px*cal distribu¬ 
tions, PYTHIA-TUNE A is re-weighted until it accurately 
follows the measured Pt'cal spectrum. The unfolding 
is carried out in two steps. First, an average correction 


is computed. The correlation < Px^jaq — P't^cal > vs 
< P'tVjal > ^ used to extract correction factors which 
are then applied to the measured jets to obtain the cor¬ 
rected transverse momenta, P'xVjor- A raw cross sec- 

, where 


tion is defined as 


N 1 ' 


P jet 

T,COR u 


,jet 

CAL 


C An jet Av jet 
^ p T,COR CAL 


IVqqr denotes the number of jets in a given PxV’OR bin, 


APt'cor the size of the bin, Ay-^ L denotes the region 
in y J cAL considered, and C is the luminosity. Second, the 
measurements are corrected for acceptance and smear¬ 
ing effects using a bin-by-bin unfolding procedure, which 
also accounts for the efficiency of the selection criteria. 

The unfolding factors, U{ P J t* cor ) = d 2 l/dp^ HAD dy^ D ’ 

u °7up T|COR uy CAL 

are extracted from Monte Carlo and applied to the mea¬ 
sured p J T et CO R distribution to obtain the final result. The 

factor [/(p J T ot co R ) increases with Px'cor an< 4 varies be¬ 
tween 1.04 at low Px cor an< 4 1-3 at very high P'xVjor- 
A detailed study of the different systematic uncertain¬ 
ties was carried out [26j. The measured jet energies were 
varied by ±2% at low p J x* and ±3% at high p^ to ac¬ 
count for the uncertainty on the absolute energy scale in 
the calorimeter |27j . This introduces an uncertainty in 
the measured cross section which varies between ± 10 % 
at low p J T ct and ^ 40 % at high p^*. A ± 8 % uncertainty 
on the jet energy resolution introduces an uncertainty 
between ±2% at low p^* and ± 8 % at high p^*. The 
unfolding procedure was repeated using HERWIG instead 
of pythia-tune A to account for the uncertainty on the 
modeling of the parton cascades and the jet fragmenta¬ 
tion into hadrons. This translates into an uncertainty 
about ±5% at low Px*. The unfolding procedure was 
also carried out using unweighted pythia-tune a, to 
estimate the residual dependence on the Px* spectrum. 
This introduces an uncertainty of ±4% above 400 GeV/c, 
which becomes negligible at lower p^*. The quoted un¬ 
certainty on e was taken into account. The effect on the 
measured cross section is less than ±3% and negligible 
for jets with pi ^ 4 above 200 GeV/c. An additional 5.8% 
uncertainty on the total luminosity is not included. 


Figure0shows the measured cross section as a function 
of Px* compared to NLO pQCD predictions. The data 
are reported in Table I. The cross section decreases by 
more than eight orders of magnitude as p^ 4 increases 
from 54 GeV/c up to about 700 GeV/c. The NLO pQCD 
predictions are computed using the JETRAD program || 
with CTEQ6.1M PDFs [ 28 ) and the renormalization and 
factorization scales (/i/j and fip) set to /io = max(p J ^i t )/2. 

Different sources of uncertainty in the theoretical pre¬ 
dictions were considered. The main contribution comes 
from the uncertainty on the PDFs and was computed 
using the Hessian method j22j. It varies from + 2 °^ at 
low p J T ct , and for pi^ about 100 GeV/c, to + 30 % 
at high p^ 1 , dominated by the limited knowledge of the 
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FIG. 1: Measured inclusive jet cross section (black dots) as 
a function of compared to NLO pQCD predictions (his¬ 
togram). The shaded band shows the total systematic uncer¬ 
tainty on the measurement. 


gluon PDF. An increase of pn and /zp from /zo to 2/zo 
reduces the theoretical predictions by 2% at low and 
8 % at high p^ 4 . Values significantly smaller than /z 0 lead 
to unstable NLO results and were not considered. 

The theoretical prediction includes a correction fac¬ 
tor, G’had (p't* ) i that approximately accounts for non- 
perturbative contributions from the underlying event and 
fragmentation into hadrons (see Table I). Chad was es¬ 
timated, using PYTHIA-TUNE A, as the ratio between the 
nominal P'x'had distribution and the one obtained by 
turning off the interactions between proton and antipro¬ 
ton remnants and the fragmentation in the Monte Carlo. 
The parton-to-hadron correction shows a strong p J ^ 4 de¬ 
pendence and decreases as p^ 4 increases from about 1.2 
at pi^ 4 of 54 GeV/c, and 1.1 for pi^ 4 about 100 GeV/c, 
to 1.02 at high p^ 4 . The uncertainty on Chad is about 
13% at low p J ^ 4 and decreases up to 1.6% at high p^ 4 , as 
determined using HERWIG instead of PYTHIA-TUNE a. 

Figure [5] shows the ratio data/theory as a function of 
p^ 4 . Good agreement is observed in the whole range in 
p^ 4 . A \ 2 test, where all the systematic uncertainties on 
the data are considered independent but fully correlated 
across different p^ 4 bins and the uncertainty on Chad 
is also included, gives a \ 2 probability of 56%. In addi¬ 
tion, FigureEJshows the ratio of pQCD predictions using 
MRST2004 j3Cl and CTEQ6.1M PDF sets, well inside 
the theoretical and experimental uncertainties. 

In summary, we have presented results on inclusive 
jet production in pp collisions at -y/s = 1.96 TeV us¬ 
ing the kx algorithm, for jets with transverse momen¬ 
tum pip 4 > 54 GeV/c and rapidity in the region 0.1 < 
|y jet | < 0.7, based on 385 pb -1 of CDF Run II data. 
The measured cross section is in agreement with NLO 
pQCD predictions after the necessary non-perturbative 



FIG. 2: Ratio Data/Theory as a function of p^ 4 . The en¬ 
closed figure expands the region p^ 4 < 298 GeV/c. The error 
bars (shaded band) show the total statistical (systematic) un¬ 
certainty on the data. A 5.8% uncertainty on the luminosity 
is not included. The solid lines indicate the PDF uncertainty 
on the theoretical prediction. The dashed line presents the 
ratio of MRST2004 and CTEQ6.1M predictions. The dotted- 
dashed line shows the ratio of predictions with 2/zq and /zq. 


D jet 

—. ± (stat.) 

dp^dyje* 

± 

(sys.) 

Chad 

± 

(stat. 

± (sys.) 

[Ge 

V/c] 



[nb/(Ge 

V/c)l 

parton —► hadron 

54 

- 62 

(14 

.6 

±0.2+|”) 

X 

10° 

1.202 

± 

0.013 

± 

0.158 

62 

- 72 

(6.53 

± 

0-041J;™) 

X 

10° 

1.154 

± 

0.003 

± 

0.113 

72 

- 83 

(2.81 

± 

0-02±“;®g) 

X 

10° 

1.134 

± 

0.005 

± 

0.094 

83 

- 96 

(1.18 

± 

o-oil2:12) 

X 

10° 

1.113 

± 

0.006 

± 

0.077 

96 

110 

(5.04 

± 

o-o4i2;5) 

X 

io- 1 

1.098 

± 

0.004 

± 

0.066 

110 

- 127 

(2.15 

± 

om+_° Q f 2 ) 

X 

10" 1 

1.079 

± 

0.005 

± 

0.047 

127 

- 146 

(8.81 

± 

o-o5iS;22) 

X 

10~ 2 

1.064 

± 

0.003 

± 

0.037 

146 

- 169 

(3.45 

± 

o-o2l2:t?) 

X 

10 -2 

1.057 

± 

0.004 

± 

0.030 

169 

- 195 

(1.28 

± 

o-oii2;i?) 

X 

10 -2 

1.047 

± 

0.003 

± 

0.023 

195 

- 224 

(4.67 

± 

o-o2l2:S) 

X 

10“ 3 

1.043 

± 

0.003 

± 

0.018 

224 

- 259 

(1.63 

± 

o.oi!°i°) 

X 

10 -3 

1.039 

± 

0.004 

± 

0.015 

259 

- 298 

(5.08 

± 

0.06lJ : °2) 

X 

10 -4 

1.034 

± 

0.003 

± 

0.010 

298 

- 344 

(1.50 

± 

°-°3l2:IS) 

X 

10- 4 

1.030 

± 

0.005 

± 

0.008 

344 

- 396 

(3.70 

± 

o-i4l2:29) 

X 

10 -5 

1.016 

± 

0.009 

± 

0.006 

396 

- 457 

(7.50 

± 

0-5511®?) 

X 

10“ 6 

1.017 

± 

0.018 

± 

0.009 

457 

- 527 

(1.31 

± 

o.22i2;t?) 

X 

10-® 

1.009 

± 

0.003 

± 

0.019 

527 

- 700 

(1.14 

± 

0.4312®?) 

X 

10- 7 

1.018 

± 

0.002 

± 

0.016 


TABLE I: Measured inclusive jet differential cross section as 
a function of Pp 4 . An additional 5.8% uncertainty on the 
luminosity is not included. The parton-to-hadron correction 
factors, Chad(Px 4 ), are applied to the pQCD predictions. 


parton-to-hadron corrections are taken into account. 
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